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Ln  j  J  ,  ABSTRACT 

An  experimental  investigation  was  made  to  determine  the  aerodynamic 
response  of  a  wind  anemometer  to  unsteady  air  loads.   The  anemometer,  a 
perforated  drag  sphere  mounted  upon  a  cylindrical  pedestal,  sensed 
dynamic  pressure  by  a  strain  gage  balance  system. 

The  anemometer  response  was  calibrated  for  static  and  dynamic 
mechanical  inputs  and  for  steady  aerodynamic  loads  prior  to  the  un- 
steady investigation.   The  unsteady  environment  was  obtained  in  an 
oscillating  flow  wind  tunnel  which^itilizes  a  rotating  shutter  valve  to 
superimpose  a  harmonically  varying  velocity  upon  a  mean  free  stream. 

For  steady  flow  conditions  the  anemometer  drag  coefficient  remained 
essentially  constant,  apparently  independent  of  Reynolds  number.   The 
anemometer  responded  well  to  the  unsteady  air  loads  at  frequencies  up 
to  20  cycles  per  second  and  mean  dynamic  pressures  to  8.85  pounds  per 
square  foot.   Unsteady  aerodynamic  response  was  qualitatively  substan- 
tiated by  a  quasi-steady  approximation  to  the  measured  unsteady  drag. 
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page  14,  line  14:   should  read  contrast  vice  constrast 

page  17,  line  26:   should  read  shown  vice  fhown. 

page  106 ,  left  hand  sketch:   should  read  steel 
rod  vice  steel  rid. 

page  109,  bottom  line:   add  a  period  at  the  end 
of  the  sentence . 

page  110,  line  3  of  the  3rd  paragraph:  the  symbol 
w  should  be  inserted  in  order  to  read  "...the 
product  wA  is . . . " . 

page  113,  the  line  under  Equation  E . 04 :   should 
read  E . o4  vice  E,04. 
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Projected  frontal  area  of  drag  sphere 
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Note:   Symbols  associated  with  a  particular  derivation  are  defined 
when  first  used  in  the  text. 
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1 . .   Introduction 

Structures  exposed  to  ground  winds  can  experience,  depending  on 
the  geometry  and/or  motion  involved,  both  steady  and  unsteady  bending, 
torsional,  and  distortional  loads.   The  effects  of  wind  loads  upon  the 
elastic  response  of  large  structures  have  been  observed  as  the  oscilla- 
tions of  tall  slender  smoke  stacks  and  antennas,  "galloping"  of  trans- 
mission lines,  warping  motion  of  suspension  bridges,  and  more  recently 
as  the  vibrations  of  space  launch  vehicles  while  on  the  launch  pad. 
These  effects  have  been  the  subjects  of  such  investigators  as  Den 
Hartog[3]  ,  Scruton  and  Walshe  flO]  ,  Scruton  [ll]  and  Buell,  et  al  [l]. 

Many  of  the  studies  conducted  on  ground  wind  loads  have  used  wind 
tunnel  tests  in  which  the  velocity  was  steady  in  magnitude  and  direction 
and  whose  profile  was  essentially  uniform  (tunnel  boundary  layer  was 
small  compared  to  model  dimensions).   In  reality  the  ground  wind  is 
rarely  constant  in  either  magnitude  or  direction  and  the  boundary  layer 
of  a  ground  wind  can  be  several  hundred  feet  thick.  A  structure's  re- 
sponse to  a  steady  uniform  flow  can  be  quite  different  from  its  response 
to  a  flow  which  varies  in  space  and  time.   Although  correct  reproduction 
of  aerodynamic  flow  properties  as  a  prerequisite  for  meaningful  dynamic 
model  tests  has  been  recognized  for  some  time,  possibly  the  first  tests 
in  which  the  dynamics  of  the  structure  and  the  turbulent  boundary  layer 
flow  were  both  reproduced  occurred  in  1964,  directed  by  A.  G.  Davenport 
f  2  J  in  cooperation  with  J.  E.  Cermak  using  the  Army  meteorological  wind 
tunnel  at  the  Colorado  State  University. 

Any  meaningful  prediction  of  wind  loads  must  be  based  on  a  know- 
ledge of  the  wind  itself,  which  in  turn  implies  accurate  measurements 
of  the  wind.   These  measurements  are  becoming  increasingly  more  important 
as  the  studies  of  wind  loading  on  large  elastic  structures  become  more 
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detailed.   Measurement  of  fluctuating,  gusty  winds  is  difficult  with 
instruments  that  are  frequency-response  limited,  and  hence  can  only 
indicate  averages  of  the  various  components.   Recently  the  National 
Aeronautics  and  Space  Administration,  Langley  Research  Center  (NASA 
Langley),  Hampton,  Virginia,  has  developed  a  fast  response  anemometer 
capable  of  following  wind  fluctuations  in  the  range  of  zero  to  five 
cycles  per  second.   According  to  Reed  and  Lynch  j[9]  this  development 
resulted  from  the  gratifying  performance  of  a  ping-pong  ball,  attached 
to  a  single  component  force  pick-up,  which  was  used  to  measure  time 
fluctuations  of  low  level  dynamic  pressures  in  the  early  1950' s.   Ad- 
ditionally, wind  tunnel  testing  at  NASA/Langley  demonstrated  that  the 
steady  drag  force  vector  for  a  perforated  hollow  sphere  remained 
aligned  with  the  steady  wind  velocity  vector  within  a  tolerance  of  +  3 
degrees.   In  constrast  a  smooth  sphere  exhibited  an  abrupt,  random 
shift  of  the  drag  vector  to  either  of  two  positions,  symmetrically 
located  about  15  degrees  from  the  wind  direction.   The  latter  result 
is  not  surprising  since  the  separation  point  for  the  flow  is  not  geo- 
metrically stabilized  as  is  the  case  for  a  perforated  sphere.   The  choice 
of  a  perforated  sphere  for  use  in  a  drag  anemometer  is  a  logical  one 
based  on  these  observations. 

The  drag  sphere  anemometer,  in  addition  to  fast  response,  has 
several  other  important  characteristics: 

(1)  measurable  drag  forces  are  generated  in  the  low  dynamic 
pressure  regions. 

(2)  the  coefficient  of  drag  is  constant  over  a  wide  range  of 
velocities,  hence  dynamic  pressure  can  be  readily  interpreted 
from  the  drag. 
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(3)  the  natural  frequency  is  high  compared  to  the  frequencies  of 
interest. 

(4)  the  drag  force  vector  remains  aligned  with  the  wind  direction. 
Current  use  of  the  drag  sphere  anemometer  has  been  indicated  by 

Duncan  and  Foughner  ["  4^|  as  part  of  a  continuing  ground  wind  loads  in- 
vestigation at  Wallops  Island,  Virginia. 

An  oscillating  flow  may  be  defined  to  consist  of  a  harmonically 
varying  velocity  superimposed  upon  a  constant  free  stream.   Low  speed 
oscillating  flow  investigations  have,  in  general,  been  limited.   Miller 
[7j  has  investigated  boundary  layer  transition  over  a  flat  plate. 
Murphy  Q8]  recently  investigated  the  pressure  distribution  about  a 
circular  cylinder  in  an  unsteady  flow.   No  investigations  concerning 
perforated  spheres  in  an  oscillating  flow  were  found  in  the  literature. 

This  investigation  was  undertaken  to  determine  the  aerodynamic 
behavior  of  the  perforated  sphere  in  an  oscillating  flow.   The  frequency 
range  of  interest  was  from  zero  to  twenty  cycles  per  second,  with  dynamic 
pressures  to  8.85  pounds  per  square  foot,  which  approximately  corresponds 
to  air  velocities  of  88  feet  per  second,  or  60  miles  per  hour. 

2.    Experimental  Equipment  and  Procedures 

The  perforated  sphere  anemometer  used  in  these  investigations  was 
supplied  to  the  Naval  Postgraduate  School  by  the  Aeroelastic  Branch  of 
NASA  (Langley)  in  support  of  the  aerodynamic  frequency  response  program. 
The  anemometer  was  floor  mounted  when  installed  in  either  the  3.5  x  5.0 
foot  low-speed  wind  tunnel  or  the  2.0  x  2.0  foot  nonsteady  wind  tunnel, 
both  of  which  are  located  at  the  Naval  Postgraduate  School,  Monterey, 
California. 
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The  seven-inch  drag  sphere  was  made  of  plastic  (weight  and  corrosion 
considerations)  and  had  forty  0.75- inch  diameter  perforations  uniformly 
distributed  over  the  surface  (see  Figure  1).   It  was  mounted  on  a  two- 
component  strain  gage  balance  system  oriented  such  that  the  balances 
sensed  the  two  horizontal  components  of  sphere  drag  force.   Static  in- 
ertial  compensation  was  provided  by  an  internal  counterweight  arrange- 
ment which  could  be  positioned  such  that  the  center  of  gravity  of  the 
sphere  system  corresponded  to  the  effective  center  of  the  strain  gage 
balances.   No  adjustments  of  the  balance  weight  were  made  following 
receipt  of  the  anemometer  from  NASA/Langley . 

Anemometer  Calibrations 

The  anemometer  strain  gage  balances,  conventional  350  ohm  full- 
bridge  circuits,  were  statically  calibrated  through  force  ranges  that 
corresponded  to  gage  strain  levels  of  approximately  +  1000  micro- inches 
per  inch.   With  a  bridge  excitation  level  of  10.0  volts  (D.C.),  the 
sensitivities  were  found  to  be  approximately  4.42  and  4.15  millivolts 
per  pound  respectively,  with  the  results  being  linear  in  the  range  in- 
vestigated.  See  Figure  2  for  actual  calibration  curves. 

Although  both  balances  were  calibrated  in  order  to  allow  simul- 
taneous evaluation  of  lateral  (side)  and  drag  forces,  only  one  balance 
was  operative  following  the  calibrations,  since  the  second  balance  be- 
came unusable  due  to  an  apparent  unbonding  of  one  of  the  strain  gage 
elements.   Consequently  the  scope  of  the  program  became  restricted  to 
a  frequency  response  evaluation  of  sphere  drag  force. 

Astrodata  differential  D.C.  amplifiers,  model  A885,  were  used  in 
these  calibrations,  with  a  gain  of  100  to  increase  the  output  level  to 
a  value  suitable  for  magnetic  tape  recording,  and  to  allow  biasing  of 
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the  steady  drag  component  so  that  only  the  unsteady  portion  of  the  signal 
would  be  recorded.   Balance  sensitivities  stated  in  a  previous  para- 
graph are  equivalent  to  amplifier  input  values. 

Since  the  anemometer  was  to  be  exposed  to  velocity  perturbations 
while  attached  to  the  floor  of  the  nonsteady  tunnel,  which  definitely  is 
a  vibratory  environment,  it  was  necessary  to  consider  the  inertial 
response  traits  of  the  anemometer  for  possible  correctional  terms 
during  data  reduction.   The  need  for  such  considerations  can  be  visual- 
ized by  realizing  that  the  strain  gage  balance  is  an  elastic  system 
that  relies  upon  relative  motion  between  the  sphere  and  the  pedestal 
base  in  order  to  produce  measurable  strains  in  the  balance  system. 

It  would  normally  be  necessary  to  consider  both  rotational  and 
translational  inertia  terms,  however,  during  preliminary  tests  the  non- 
steady  wind  tunnel  (described  later  in  this  section)  was  found  to  have 
floor  motion  predominantly  in  a  direction  tangential  to  the  floor  with 
maximum  accelerations  of  2.0  "g's"  at  10.0  and  17.5  cps.   This  accelera- 
tion direction  is  dominant  since  the  tunnel  test  section  is  relatively 
stiff,  and  the  tunnel  is  mounted  on  a  framework  of  steel  pipes  which 
act  like  a  "bent",  without  significant  bending  moment  constraint  at 
the  joints.   Therefore  a  calibration  was  made  only  in  the  translational 
vibration  mode. 

The  inertial  frequency  response  calibration  was  performed  over  a 
frequency  range  of  5  to  100  cps.   A  Calidyne  Shaker,  model  219,  operated 
by  a  Ling  control  console,  model  CP  3/4  was  used  to  accelerate  a  steel 
plate  upon  which  the  anemometer  and  reference  accelerometer  were  mounted, 
as  fhown  in  Figure  1.   Horizontal  accelerations  were  measured  with  a 
Statham  accelerometer,  model  SCO  210,  which  was  used  as  a  reference 
sensor  unit  in  both  the  mechanics  laboratory  and  the  nonsteady  wind 
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tunnel.   The  primary  range  of  frequencies  was  below  100  cps,  and  this 
was  readily  measured  by  the  Statham  strain  gage  device  since  it  had 
a  usable  frequency  range  from  0(D.C.)  to  150  cps,  and  an  acceleration 
range  of  +  5  "g's"  (see  Appendix  A  for  accelerometer  calibration).   The 
shaker  was  connected  to  the  plate  by  a  1/16  inch  diameter  steel  wire 
flexure  coupling,  which  had  an  estimated  critical  buckling  load  of 
5400  pounds.   A  wire  flexure  coupling  was  used  to  avoid  imposing  trans- 
verse loads  onto  the  shaker  table.   The  steel  plate  rested  upon  a  steel 
surface  table  which  was  adjusted  to  be  level.   The  surface  table  and 
the  plate  had  polished  surfaces,  and  an  oil  film  was  applied  between 
the  moving  surfaces  in  order  to  reduce  the  level  of  frictional  forces. 
The  coupling  was  positioned  such  that  the  axial  force  transmitted 
through  the  wire  coupling  acted  approximately  through  the  center  of 
mass  of  the  combined  dynamic  assembly  (plate,  anemometer,  and  accelero- 
meter) . 

If  the  anemometer  were  represented  as  a  system  with  a  single  degree 
of  freedom:  n i  • 


1 


Li) 


then,  the  equation  of  motion  for  the  freely  vibrating  mass  is 
Let  H  =  X-  W   ,  and  rewrite 

bo  &   4-  CH  ■+  ki   =  -"^3 


then  for  sinusoidal  motion  of  the  base 


Stn  Ojt 


(2.01) 


(2.02) 


18 


where  , 


2^  =  -^-    j 


The   solution  of    (2.02)    is: 

[i  -  sr] 

where  2~  (t)  is  the  displacement  of  the  anemometer  force  pickup  rela- 
tive  to  the  base,  and  M«>uJ  is  the  magnitude  of  the  harmonic  forcing 
function. 

Accelerations  of  various  magnitudes  and  constant  frequency  were 
applied  to  the  base  plate  at  four  intervals  up  to  and  down  from  the 
maximum  of  2.0  "g's".   The  anemometer  response  £  (t) ,  accelerometer 
output,  W0U*  ,  and  phase  difference  between  these  sensors  were  recorded 
at  each  level  of  acceleration.   The  frequency  was  increased  and  the 
procedure  repeated.   Figure  3  is  a  summary  of  the  anemometer  inertial 
response  calibration. 

The  anemometer  natural  frequency  was  found  to  be  approximately 
60  cps  as  indicated  by  maximum  response  from  oscilloscope  traces,  and 
by  observation  of  large  excursions  of  the  perforated  sphere.   For 
frequencies  greater  than  60  cps,  the  response  wave  forms  became  very 
distorted.   Considerable  difficulty  was  encountered  in  the  determination 
of  a  90  degree  phase  shift  between  input  and  output.   The  difficulty 
was  attributed  to  the  anemometer  not  behaving  as  a  system  with  a  single 
degree  of  freedom. 
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Anemometer  damping  was  evaluated  by  observing  the  output  response 
to  a  step  input.   No  low  frequency  oscillation  ensued  for  any  of  the 
step  inputs,  and  it  was  concluded  that  the  system  fundamental  mode  was 
at  least  critically  damped.   Examination  of  the  response,  however,  with 
the  single  sweep  capability  of  the  Tektronix  (type  551)  oscilloscope, 
revealed  a  high  frequency  (greater  than  100  cps)  signal  superimposed 
upon  the  basic  wave  form.   The  excitation  of  this  high  frequency  signal 
was  controllable  to  a  degree  and  was  apparently  caused  by  a  ringing 
mode  of  the  perforated  shell. 

Nonsteady  Tunnel  Investigation 

The  investigation  of  anemometer  unsteady  aerodynamic  response  was 
conducted  in  the  low  speed,  oscillating  flow  tunnel  located  in  the 
Aeronautical  Laboratories  of  the  Naval  Postgraduate  School.   The  open 
circuit  tunnel  has  a  24  inch  square  test  section,  223.5  inches  long. 
Air  is  driven  by  two  series-connected  Joy  Axivane  fans  connected 
directly  to  100  hp,  1750  rpm  electric  motors.   The  tunnel  has  a  bell- 
mouthed  entrance  with  a  nozzle  area  contraction  ratio  of  16:1.   There 
are  three  fine  mesh  screens  in  the  settling  section  area  which  reduce 
steady  flow  turbulence  levels  in  the  test  section  below  0.1  percent 
of  mean  velocity. 

Construction  of  the  test  section  floor  and  ceiling  is  from  2.0  inch 
thick  aluminum,  the  sides  from  2.0  inch  lucite.   Heavy  construction  was 
necessary  to  reduce  structural  deflections  caused  by  large  variations 
in  static  pressure. 

Oscillating  flow  was  created  by  a  rotating  shutter  valve  assembly 
located  at  the  trailing  edge  of  the  test  section  and  driven  through 
belts  by  a  variable  speed  electric  motor.   The  shutter  blades  were 
removable  to  permit  installation  of  various  widths  that  could  enclose 
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from  50  to  100  percent  of  the  test  section  area,  and  which  provided 
unsteady  velocity  amplitudes  from  1.3  to  84.0  percent  of  free  stream 
velocity.   Various  pulley  ratios  of  the  belt  drive  were  available 
effecting  an  operational  range  of  frequency  from  2  cps  to  933  cps. 

Murphy  £  8j  performed  steady  flow  velocity  surveys  across  the  test 
section  and  found  that  the  variation  of  V    relative  to  Vmax  was  less 
than  +  1  percent  for  distances  greater  than  three  inches  from  the 
tunnel  walls. 

For  the  actual  nonsteady  tunnel  tests  it  was  desired  to  subject 
the  aneometer  to  five  values  of  mean  dynamic  pressure  with  each  of  the 
three  available  shutter  vane  sizes.   The  investigation  in  general  in- 
cluded fifteen  points  of  data  for  each  run  at  various  frequencies  up 
to  20  cps. 

In  order  to  determine  the  effect  of  blockage  produced  by  the  ro- 
tating shutter  vanes,  and  to  insure  that  the  hot  wire  system  would  per- 
form satisfactorily,  several  trial  runs  were  performed  with  the  anemom- 
eter absent  from  the  test  section.   Dynamic  pressures  of  7.8  and  15.6 
pounds  per  square  foot  (psf,)  were  set  with  the  vanes  horizontal.  The 
hot  wire  DC  output  was  set  to  read  full  scale  (1.0  volts)  for  zero 
frequency  and  maximum  velocity.   The  tunnel  was  operated  at  various 
frequencies  and  the  drop  of  the  DC  level  of  the  hot  wire  was  noted. 
Figure  4  is  typical  of  results  for  an  initial  mean  dynamic  pressure  of 
7.8  psf.   Post-run  checks  were  performed  to  ensure  the  absence  of  hot 
wire  drift. 

With  the  information  obtained  from  an  evaluation  of  the  blockage 
effect,  a  zero  frequency  dynamic  pressure  was  selected.  The  hot  wire 
DC  attenuation  was  adjusted  to  full  scale,  and  the  zero  frequency  RMS 
values  were  recorded  for  the  anemometer,  hot  wire,  and  accelerometer. 
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The  speed  of  the  rotating  shutter  valve  was  initially  adjusted  to  the 
minimum,  3.6  cps,  for  the  combination  of  pulley  ratios  used.   After 
amplified  anemometer  and  hot  wire  DC  values  were  biased  to  zero  and 
the  oscillation  frequency  was  stabilized,  tape  recording  of  signals 
for  the  particular  data  point  was  commenced.   RMS  values  of  the  ane- 
mometer, hot  wire,  frequency,  and  mean  dynamic  pressure  were  logged 
manually.   This  procedure  was  repeated  at  seven  increments  of  frequency 
up  to  and  down  from  the  maximum  of  20  cps.   Post-run  zero  frequency  and 
zero  dynamic  pressure  were  noted.   This  procedure  of  data  acquisition 
was  repeated  at  five  levels  of  dynamic  pressure  for  each  of  the  shutter 
vane  widths. 

Nonsteady  Tunnel  Instrumentation 

The  Ampex  magnetic  tape  system  used  was  an  electron  tube  type 
model  FR-100  with  14-channel  recording  capability.   Both  frequency 
modulated  and  direct  record  amplifiers  were  available.   A  tape  speed 
of  7.5  inches  per  second  was  used  throughout  the  data  acquisition, 
which  provided  a  flat  frequency  response  on  the  FM  channels  from  0  to 
1250  cps.   Three  frequency  modulated  channels  were  used  to  record  ane- 
mometer, hot  wire,  and  accelerometer  signals,  whereas  voice  information, 
to  identify  blocks  of  taped  data,  was  recorded  through  a  direct  record 
amplifier. 

A  Tektronix  type  555  oscilloscope  with  a  pair  of  4-channel  pre- 
amplifiers was  used  to  simultaneously  monitor  original  and  recorded 
data. 

An  integrating  digital  voltmeter  (Vidar  510)  was  used  to  alter- 
nately monitor  the  DC  output  of  the  anemometer  bridge  circuit  and  the 
hot  wire  anemometer.   The  DC  values  of  the  amplified  signals  from  both 
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instruments  were  maintained  at  values  as  near  to  zero  as  possible  by  a 
battery-supplied  bias  voltage.   The  DC  biasing  permitted  the  tape  re- 
cording of  optimum  level  AC  signals.   Figure  5  presents  a  block,  diagram 
of  the  entire  data  acquisition  system  shown  by  photograph  in  Figure  6. 
Details  of  the  anemometer  installation  are  shown  in  Figure  7.   Shielded 
co-axial  cables  with  BNC  fittings  were  used  to  reduce  the  influence  of 
extraneous  signals,  for  example  60-cycle  ground  loops.   A  centrally 
located  switch  panel  afforded  an  orderly  flow  of  data  and  convenient 
monitoring  of  information. 

The  hot  wire  probe  was  constructed  at  the  Postgraduate  School  by 
soldering  an  0.0015  inch  diameter  tungsten  wire  on  a  pair  of  jewelers 
broaches.   The  jewelers  broaches  were  imbedded  in  epoxy  plastic,  which 
was  fitted  to  a  1/4  inch  hollow  steel  rod.   The  performance  of  the 
Security  Associates  bridge  circuit/analog  computer  system  was  excellent. 
Repeated  results  assured  a  linear  DC  variation  from  zero  to  full  scale 
and  a  high  degree  of  system  stability  was  noted.   The  sensitivity  of 
the  system  was  such  that  outside  wind  fluctuations  during  steady  opera- 
tion of  the  wind  tunnel  were  detectable. 
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3.   Experimental  Results  and  Discussion: 
Run  Program 
The  investigation  of  the  aerodynamic  response  of  the  wind  anemometer 
performed  in  the  oscillating  flow  tunnel  was  separated  into  three  phases: 
(a,)  tunnel  familiarization  and  calibration  of  flow,  (b.)  evaluation 
of  the  steady  flow  drag,  and  (c.)  the  determination  of  unsteady  drag 
in  oscillating  flow.   The  phases  may  be  further  subdivided  as  follows: 
Phase  I    Calibration 
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Data  Reduction 

The  initial  step  in  the  reduction  of  experimental  data  was  to 
evaluate  the  operating  conditions  for  which  the  anemometer  output  was 
significantly  affected  by  tunnel  floor  accelerations.   In  consideration 
of  the  inherent  fluctuations  of  anemometer  RMS  output,  "significantly" 
was  interpreted  to  mean  those  data  points  which  influenced  the  total 
anemometer  signal  by  more  than  five  percent.   Frequencies  which  caused 
excessive  tunnel  vibration  were  avoided,  but  using  the  criteria  above 
it  was  necessary  to  correct  about  25  percent  of  the  data  points. 

A  method  to  correct  the  anemometer  response,  both  in  magnitude 
and  phase,  is  presented  in  Appendix  B.  An  important  assumption  in  the 
analysis  is  that  the  dominant  floor  accelerations  are  of  a  sinusoidal 
nature  and  of  the  same  frequency  as  the  oscillating  air  flow.   The 
assumption  was  based  upon  oscilloscope  photographs  that  were  made  of 
the  response  signals  reproduced  from  the  tape  recorded  data.   Two  Kron- 
Hite  band  pass  filters  were  employed  in  the  playback  channels  to  elim- 
inate the  high  frequency  components  of  the  actual  signals.   Tables  1, 
2  and  3  contain  the  corrected  data  for  small,  medium  and  large  shutter 
vanes  respectively. 

Appendix  C  is  a  development  of  hot  wire  response  to  a  base  or 
floor  acceleration  based  upon  elementary  beam  theory.   For  all  values 
of  frequency  and  acceleration  in  the  range  investigated  this  effect 
can  be  neglected.   The  effect  of  tunnel  floor  velocity  on  hot  wire 
response,  considered  in  Appendix  D,  is  less  than  two  percent  of  the 
unsteady  air  velocity  and  therefore  negligible.   The  hot  wire  output 
can  then  be  considered  as  the  true  behavior  of  the  oscillating  air 
flow. 
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The  phase  differences  between  (a.)  anemometer  and  hot  wire  signals 
and  (b.)  anemometer  and  accelerometer  signals  were  determined  by 
averaging  the  output  of  an  AD-YU  phase  meter  for  an  interval  of  taped 
data  (30-45  sees).   The  difference  was  qualitatively  substantiated  by 
photographic  measurements.   Phase  angles  of  the  anemometer  relative  to 
the  hot  wire  are  listed  in  Tables  4  and  5  for  the  small  and  medium 
shutter  vanes.   Phase  differences  for  the  large  vanes  are  omitted  due 
to  the  large  degree  of  variation  resulting  from  the  complex  waveform. 

Drag-Steady  Values 

Figures  8  and  9  are  plots  of  raw  anemometer  drag  as  a  function  of 
dynamic  pressure  as  obtained  in  the  steady  West  Coast  Research  (W.C.R.) 
tunnel  and  in  the  oscillating  flow  tunnel  for  zero  frequency.   The 
steady  coefficient  of  drag  was  evaluated  for  each  set  of  data  using 
the  anemometer  static  sensitivity  (Figure  2)  previously  determined. 
The  drag  coefficient  obtained  from  the  W.C.R.  tunnel  was  0.69,  while 
the  non-steady  tunnel  yielded  0.71,  which  includes  an  estimated  cor- 
rection for  solid  blocking.   For  both  evaluations,  the  anemometer  out- 
put was  a  linear  function  of  dynamic  pressure  in  the  range  investigated, 
which  indicated  that  the  drag  coefficient  for  the  perforated  drag  sphere 
was  apparently  independent  of  Reynolds  Number. 

Reed  L.  9  J  has  shown  for  the  one-dimensional  case  of  a  non-steady 
velocity  v(t)  superimposed  upon  a  mean  free  stream  V  that  the  time 
average  of  the  instantaneous  velocity  squared  is 


vV/O+JJc) 
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Taking  the  square  root  of  both  sides,  expanding  the  right  hand 
side  in  a  binomial  series,  and  retaining  first  order  terms  yields: 

In  other  words,  the  mean  velocity  when  determined  by  a  measurement  of 

2 

dynamic  pressure  (or  V  )  is  biased  by  the  turbulence  in  the  flow. 

The  unsteady  tunnel  was  used  to  determine  the  effect  of  an  oscil- 
lating flow  on  the  mean  drag  of  the  anemometer.   The  first  step  was  to 
measure  the  drag  in  a  steady  flow  over  a  dynamic  pressure  range  of 
1.04  to  15.6  pounds  per  square  foot.   This  test  corroborated  earlier 
results  obtained  in  the  W.C.R.  tunnel.   Next  the  mean  anemometer  drag 
was  recorded  in  an  oscillating  flow  (both  small  and  medium  sized  vanes 
were  used)  over  a  mean  dynamic  pressure  range  of  2.08  to  8.84  pounds 
per  square  foot.   The  mean  drag  was  recorded  at  about  16  discrete 
frequencies  between  3.6  and  19.1  cycles  per  second  at  each  of  9  dynamic 
pressure  levels.   It  was  found  that  the  mean  anemometer  drag  at  a  mean 
dynamic  pressure  was  independent  of  flow  oscillation  frequency  and 
equal  to  the  steady  flow  value.   The  results  of  the  steady  and  the  un- 
steady flow  tests  are  shown  in  Figure  9. 

Unsteady  Drag  and  Phasing 

The  variation  of  unsteady  drag  as  a  function  of  frequency  and 
dynamic  pressure  for  the  case  of  the  small  shutter  vanes  is  presented 
in  Figure  10.   The  individual  data  points  for  the  curves  were  obtained 
by  positive  incremental  frequency  steps  up  to  the  maximum  obtainable, 
and  then  by  negative  steps  to  the  minimum  frequency.   A  high  degree  of 
repeatability  was  obtained  for  the  small  vanes.   Unsteady  drag,  for 
the  lowest  dynamic  pressure  (1.0  psf)  in  the  series,  remained  essentially 
constant  with  frequency,  while  each  of  the  other  curves  for  higher 
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dynamic  pressures  generally  decreased  with  frequency  and  apparently 
approached  a  limiting  value.   It  is  interesting  to  note  that  the  varia- 
tion of  unsteady  drag  with  dynamic  pressure  at  a  constant  frequency 
of  19  cps  is  approximately  one-fifth  the  value  at  a  frequency  of 
3.6  cps. 

All  of  the  curves  exhibited  a  peak  in  the  vicinity  of  10  cps, 
the  relative  magnitude  of  which  increased  with  dynamic  pressure.   Also, 
at  frequencies  near  10  cps,  the  ratio  of  v/V  was  a  local  maximum  for 
the  particular  run,  and  is  believed  to  cause  the  rise  in  unsteady  drag. 
The  reason  for  an  increase  in  this  velocity  perturbation  ratio  at  this 
frequency  is  not  known  but  was  attributed  to  the  flow  characteristics 
of  the  tunnel. 

The  unsteady  drag  variation  for  the  medium  and  large  shutter 
vanes  is  shown  in  Figures  11  and  12.   For  the  case  of  medium  vanes 
an  additional  peak  appears  at  approximately  17.5  cps,  which  also  can 
be  attributed  to  a  local  maximum  of  v/V.   Although  the  medium-vane 
curves  display  a  lower  degree  of  consistency  than  do  the  small-vane 
curves,  the  trend  of  decreasing  magnitudes  with  frequency  is  maintained, 
as  is  the  trend  toward  a  limiting  value  for  each  dynamic  pressure  with 
increasing  frequency.   The  large  shutter-vane  results  display  an  even 
lower  degree  of  consistency,  but  still  display  the  drag  peaks  at 
approximately  10  and  17  cps.   Unsteady  drag  dependence  on  dynamic 
pressure,  as  indicated  by  the  smaller  vane  data,  is  substantiated  for 
frequencies  less  than  11  cps;  for  frequencies  greater  than  11  cps  the 
scatter  of  data  precludes  immediate  interpretation. 

The  anemometer  phase  lag  to  an  aerodynamic  forcing  function  is 
affected  both  by  dynamic  pressure  and  frequency  as  shown  in  Figure  13, 
which  is  a  summary  of  results  for  the  small  velocity  perturbations 
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(v/V  =  0.02).   Phase  lag  generally  increases  with  frequency  and 
apparently  has  no  dependence  upon  dynamic  pressure.   The  increase  in 
phase  lag  with  frequency  up  to  about  16  cps  was  also  observed  during 
the  inertial  investigation  described  in  Section  2  (cfi,  Figure  3). 

With  the  exception  of  the  lowest  dynamic  pressure  (0.73  psf)  in 
the  medium  vane  runs  (where  v/V  =  0.12),  anemometer  phase  lag  re- 
mained near  zero  up  to  a  frequency  of  12  cps,  where  a  rapid  increase 
was  observed  toward  a  peak  value  occurring  at  approximately  15  cps. 
See  Figure  14. 

Quasi- Steady  Assumption 

A  theoretical  treatment  of  bluff  body  oscillating  flows  is  very 
scarce  in  published  literature.  Isaacs  £ 6  J  ,  in  an  investigation  of 
the  non-steady  flow  over  a  helicopter  blade  at  small  angle  of  attack, 
has  approached  the  problem  for  a  potential  flow  case.   However,  the 
analysis  included  two  limiting  assumptions:   the  flow  was  two-dimen- 
sional over  a  flat  plate  airfoil,  and  the  vortices  shed  remain  in  the 
plane  of  the  airfoil.   The  method  is  not  directly  applicable  to  the 
three-dimensional  sphere  problem  because  of  obvious  differences  be- 
tween idealized  potential  flow  and  bluff  body  separated  flow,  but 
Isaacs  does  suggest  that  the  lift  on  the  airfoil  can  be  evaluated 
from  an  averaging  of  the  instantaneous  lift  due  to  the  variable  velocity, 

Fung  [^ 5]  further  defines  this  quasi-steady  principle,  which  assumes 
that  the  characteristics  of  an  airfoil  in  an  unsteady  flow  at  any  in- 
stant of  time  are  equal  to  the  characteristics  of  the  same  airfoil  in 
a  steady  flow  with  flow  parameters  equal  to  the  instantaneous  values. 

Appendix  E  is  a  development  of  the  quasi-steady  principle  for  the 
three-dimensional  sphere. 
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Drag  Ratio 
Anemometer  drag  ratio  (DR^)  is  defined  as  the  quotient  of  un- 
steady drag  divided  by  drag  produced  by  the  steady  component  of  velocity, 

A 

The  quasi-steady  drag  ratio  (DRq)  is  approximately  equal  to  2  —  and  is 

developed  in  Appendix  E.   These  ratios  are  plotted  versus  Strouhal 

f   i 
number  (-=—=■)  to  indicate  the  variation  with  frequency  and  mean  velocity, 

and  to  evaluate  the  degree  of  validity  of  the  quasi-steady  assumption. 

The  scales  for  the  coordinate  axes  are  not  consistent  due  to  the  large 

variation  of  maximum  values  between  runs. 

The  quasi-steady  assumption  agrees  best  with  the  anemometer 

response  for  the  runs  with  the  small  shutter  vanes  installed.   Two 

factors  are  significant  for  this  result:   (1)  velocity  perturbations 

produced  by  the  small  vanes  are  better  approximations  to  harmonic 

A  k 

functions,  and  (2)  the  ratio  of  v/V  is  small,  which  is  inherent  to 
the  approximation  of  quasi-steady  theory.   Photographs  of  response 
signals  from  oscilloscope  traces  are  presented  in  Figures  15  and  16  to 
illustrate  the  effects  of  shutter  vane  size  upon  the  deviation  from 
harmonic  perturbations. 

The  variation  of  drag  ratio  with  Strouhal  number  and  dynamic 
pressure  for  five  runs  with  the  small  shutter  vanes  installed  is  con- 
tained in  a  sequence  of  figures  beginning  with  Figure  18.   The  ane- 
mometer response  curves  have  the  same  form  as  the  curves  for  the  quasi- 
steady  approach  (based  upon  the  hot  wire  output)  or,  in  other  words, 
the  anemometer  is  sensing  the  velocity  perturbations  throughout  the 
frequency  range.   The  level  of  response  increases,  as  expected,  with 
increasing  dynamic  pressure.   However,  for  a  given  dynamic  pressure, 


30 


both  the  average  value  of  the  drag  ratio  and  the  quotient,  DR./DRq, 
decrease  with  increasing  Strouhal  number.   Figures  20,  21  and  22  also 
contain  results  obtained  from  separate  test  dates  to  demonstrate  the 
repeatability  of  data  and  thereby  provide  a  confirmation  of  results. 

Drag  ratio  variation  obtained  with  the  medium  vanes  installed  is 
presented  in  a  second  sequence  of  figures,  starting  with  Figure  23. 
The  correlation  of  anemometer  and  quasi-steady  results  is  still  main- 
tained, although  the  appearance  of  a  second  peak  occurs  at  a  Strouhal 
number  of  0.22  (the  equivalent  frequency  is  approximately  17.5  cps). 
For  the  medium  vanes,  the  average  value  of  the  drag  ratios  are  on  the 
order  of  seven  times  the  small  vane  values,  although  the  actual  vane 
size  has  been  increased  only  50  percent. 

Large  vane  effects  upon  drag  ratio  exhibit  more  deviation  and 
less  repeatability  than  the  vane  sizes  previously  considered.   The 
large  vanes  create  more  distortion,  of  the  velocity  perturbation,  and 
anemometer  response  becomes  more  random  in  nature.   Large  fluctuations 
(50%  of  the  freestream  value)  produced  anemometer  response  that  was 
generally  about  half  of  the  quasi-steady  value,  as  shown  in  the  series 
of  figures  starting  with  Figure  28. 

Due  to  the  non-linear  variation  of  the  tunnel  nonsteady  velocity 
for  a  given  mean  dynamic  pressure  and  shutter  size,  as  a  function  of 
frequency,  it  is  not  clear  what  the  effects  are,  if  any,  of  aerodynamic 
damping  upon  anemometer  response.   However,  a  tendency  toward  reduced 
response  (as  evidenced  by  the  reduced  drag  ratio  in  relation  to  the  hot 
wire  value)  has  been  observed  for  a  combination  of  the  greater  dynamic 
pressures,  large  fluctuations,  and  higher  frequencies. 
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Throughout  the  nonsteady  investigation,  the  tunnel  radiated 
various  levels  of  noise  dependent  upon  frequency,  shutter  size,  and 
dynamic  pressure.   If  the  inlet  test  section  is  compared  to  an  open 
organ  pipe,  then  the  fundamental  frequency  may  be  calculated  equal  to 
27  cycles  per  second.   The  peaks  of  the  drag  ratio  curves  may  then  re- 
flect a  significant  degree  of  sub-harmonic  response. 

Harmonic  Considerations 

The  output  signals  obtained  from  the  anemometer  and  hot  wire  were 
obviously  not  true  sinusoids.   Observations  of  these  signals  indicated 
a  possible  significant  higher  harmonic  component,  as  shown  in  Figure  33. 
The  anemometer  response  has  been  shown  to  follow  the  gross  aerodynamic 
input  reasonably  well  but  a  more  quantitative  definition  of  response 
was  desired.   Therefore,  the  gross  signals  were  examined  to  determine 
the  amplitude  and  phase  relationships  among  the  fundamental,  second, 
and  third  harmonic  components  and  between  the  anemometer  and  hot  wire 
signals.   A  numerical  approximation  of  a  Fourier  analysis  was  used,  as 
indicated  in  Appendix  F.   Since  the  required  computations  were  performed 
with  a  desk  calculator,  only  a  few  selected  data  points  were  analyzed. 
Results  of  this  examination  are  shown  in  Table  6.   As  a  check  on  the 
consistency  of  the  approximations,  and  the  results,  two  different  data 
samples  from  runs  1-2  and  2-4  were  analyzed.   The  data  samples  from 
run  1-2  were  taken  from  different  cycles  on  the  same  oscilloscope 
photograph,  while  the  samples  from  run  2-4  were  taken  from  separate 
photographs.   Quantities  with  good  correlation,  and  discrepancies,  are 
obvious  in  Table  6.   The  existence  of  a  significant  second  harmonic 
component  is  definitely  shown,  having  an  average  amplitude,  over  the 
points  tested,  of  23  percent  and  17  percent  of  the  fundamental  component 
amplitudes  for  the  hot  wire  and  anemometer  respectively.   Corresponding 
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values  for  the  third  harmonic  component  are  4  and  6  percent.   The 
fundamental  component  of  the  hot  wire  led  that  of  the  anemometer  by 
an  average  of  9  degrees.   The  phase  angles  between  hot  wire  and  ane- 
mometer computed  for  the  second  and  third  harmonic  components  were 
random  in  sign  and  magnitude. 

Inspection  of  several  oscilloscope  photographs  (cf.  Figure  34) 
revealed  a  higher  frequency  superimposed  on  the  basic  signal.   This 
higher  frequency,  about  145  cps,  shows  up  on  both  hot  wire  and  ane- 
mometer signals,  but  it  is  much  more  prominent  for  the  latter.   The 
sensitivity  of  the  anemometer  to  this  input  frequency  is  believed  due 
to  the  oscillations  of  the  sphere  itself,  as  opposed  to  the  sphere- 
pedestal  combination.   This  "ringing"  of  the  sphere  was  also  noted 
during  part  of  the  inertial  response  calibration  to  determine  ane- 
mometer damping  ratio. 
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4.   Conclusions 

The  natural  frequency  of  the  anemometer  was  60  cps  with  approxi- 
mately a  critical  amount  of  damping. 

For  steady  flow  conditions  the  anemometer  has  a  constant  drag 
coefficient,  thus  apparently  independent  of  Reynolds  Number  within 
the  dynamic  pressure  range  of  0.78  to  15.6  psf. 

The  average  drag  measured  by  the  anemometer  in  an  unsteady  flow 
for  a  given  mean  dynamic  pressure  was  independent  of  frequency  and 
equal  to  the  drag  measured  in  a  steady  flow  at  the  same  dynamic 
pressure. 

Unsteady  drag,  as  determined  from  a  quasi-steady  approach,  pro- 
vided a  good  qualitative  comparison  to  the  unsteady  drag  measured 
by  the  anemometer.   Quantitative  results  of  the  quasi-steady  assumption 
provided  only  an  order  of  magnitude  indication  of  anemometer  response. 

The  unsteady  component  of  anemometer  response  increased  with  mean 
dynamic  pressure  and  with  unsteady  velocity.   For  the  runs  performed 
with  the  small  shutter  vanes,  the  unsteady  anemometer  response  rela- 
tive to  the  unsteady  velocity  decreased  with  increasing  frequency, 
and  apparently  approached  a  limiting  value. 

The  investigation  was  complicated  by  an  amplification  of  the  un- 
steady air  flow  at  several  frequencies  in  the  range  of  interest,  and 
by  a  significant  second  harmonic  component,  about  25  percent  of  the 
fundamental  amplitude. 

It  is  recommended  that  future  investigation  be  directed  toward  a 
thorough  definition  of  the  oscillating  flow  behavior  of  the  nonsteady 
tunnel  throughout  its  operating  range. 
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TABLE  1 

REDUCED  DATA:  SMALL 

VANES 

A 

Run* 

f 

q 

eA 

2  X 
V 

(CPS) 

(psf) 

(VRMS) 
x  102 

1-1 

3.6 

4.21 

0.031 

0.068 

1-2 

5.0 

4.21 

0.026 

0.054 

1-3 

7.7 

4.21 

0.018 

0.042 

1-4 

10.0 

4.21 

0.022 

0.056 

1-5 

13.0 

4.21 

0.013 

0.039 

1-6 

16.3 

4.21 

0.011 

0.036 

1-7 

19.2 

4.21 

0.010 

0.044 

1-8 

17.8 

4.21 

0.011 

0.039 

1-9 

15.2 

4.16 

0.012 

0.036 

1-10 

11.6 

4.16 

0.016 

0.052 

1-11 

9.2 

4.16 

0.019 

0.048 

1-12 

6.4 

4.16 

0.020 

0.046 

1-13 

4.2 

4.16 

0.029 

0.062 

1-14 

3.6 

4.16 

0.031 

0.070 

*The  run  number  is  followed  by  a  coordination  number  which  is  used 
to  identify  data  as  recorded  on  magnetic  tape. 
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TABLE   1    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 
x  102 

2^ 
V 

2-1 

3.6 

2.34 

0.0130 

0.044 

2-2 

6.0 

2.34 

0.0100 

0.032 

2-3 

8.8 

2.34 

0.0084 

0.031 

2-4 

11.2 

2.34 

0.0090 

0.034 

2-5 

13.3 

2.34 

0.0068 

0.026 

2-6 

15.0 

2.34 

0.0068 

0.025 

2-7 

16.9 

2.34 

0.0064 

0.075 

2-8 

19.1 

2.34 

0.0064 

0.029 

2-9 

18.1 

2.34 

0.0060 

0.027 

2-10 

15.9 

2.34 

0.0065 

0.025 

2-11 

14.1 

2.34 

0.0066 

0.025 

2-12 

10.1 

2.34 

0.0099 

0.039 

2-13 

7.0 

2.34 

0.0090 

0.028 

2-14 

5.0 

2.34 

0.0110 

0.034 

2-15 

4.2 

2.34 

0.0120 

0.038 
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TABLE  1    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 
(VRMS) 

x  102 

A. 

2  % 

V 

0.077 

3-1 

3.6 

7.60 

0.087 

3-2 

5.8 

7.60 

0.050 

0.064 

3-3 

8.6 

7.55 

0.037 

0.056 

3-4 

11.1 

7.55 

0.047 

0.070 

3-5 

13.6 

7.55 

0.026 

0.054 

3-6 

17.4 

7.55 

0.017 

0.058 

3-7 

19.2 

7.60 

0.018 

0.078 

3-8 

15.7 

7.60 

0.019 

0.054 

3-9 

12.4 

7.55 

0.033 

0.058 

3-10 

9.7 

7.55 

0.050 

0.076 

3-11 

7.0 

7.55 

0.043 

0.056 

3-12 

4.8 

7.60 

0.060 

0.078 

3-13 

4.2 

7.55 

0.070 

0.078 
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TABLE   1    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 

x  102 

2^ 

V 

4-1 

3.6 

5.98 

0.055 

0.0594 

4-2 

5.1 

5.98 

0.043 

0.0458 

4-3 

7.9 

5.92 

0.030 

0.0344 

4-4 

10.1 

5.92 

0.040 

0.0504 

4-5 

12.0 

5.92 

0.030 

0.0412 

4-6 

14.0 

5.92 

0.019 

0.0322 

4-7 

17.0 

5.88 

0.016 

0.0300 

4-8 

19.3 

5.98 

0.016 

0.0346 

4-9 

15.6 

5.92 

0.017 

0.0298 

4-10 

13.1 

5.92 

0.021 

0.0344 

4-11 

9.1 

5.92 

0.034 

0.0392 

4-12 

6.6 

5.92 

0.033 

0.0390 

4-12 

4.6 

5.98 

0.046 

0.0504 

4-14 

4.0 

5.98 

0.050 

0.0574 
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TABLE  1    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 

x  102 

2^ 
V 

10-1 

3.6 

0.99 

0.0039 

0.0288 

10-2 

5.5 

0.99 

0.0037 

0.0266 

10-3 

7.4 

0.99 

0.0038 

0.0266 

10-4 

9.7 

0.99 

0.0041 

0.0308 

10-5 

12.3 

0.99 

0.0039 

0.0236 

10-6 

14.5 

0.99 

0.0030 

0.0206 

10-7 

16.7 

0.99 

0.0035 

0.0194 

10-8 

19.2 

0.99 

0.0032 

0.0184 

10-9 

17.6 

0.99 

0.0032 

0.0190 

10-10 

13.3 

0.99 

0.0036 

0.0224 

10-11 

10.9 

0.99 

0.0039 

0.0288 

10-12 

8.4 

0.99 

0.0038 

0.0266 

10-13 

6.3 

0.99 

0.0037 

0.0266 

10-14 

4.3 

0.99 

0.0039 

0.0288 
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TABLE 
REDUCED  DATA: 

2 

MEDIUM  VANES 

Run* 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 

x  102 

2  2 

V 

11-1 

3.7 

1.56 

0.048 

0.282 

11-2 

6.4 

1.45 

0.033 

0.228 

11-3 

8.6 

1.40 

0.037 

0.216 

11-4 

11.2 

1.45 

0.036 

0.228 

11-5 

13.1 

1.40 

0.032 

0.179 

11-6 

15.5 

1.35 

0.033 

0.202 

11-7 

18.2 

1.40 

0.033 

0.190 

11-8 

19.3 

1.51 

0.029 

0.166 

11-9 

16.8 

1.35 

0.033 

0.212 

11-10 

14.2 

1.35 

0.033 

0.184 

11-11 

9.7 

1.40 

0.040 

0.278 

11-12 

717 

1.40 

0.030 

0.222 

11-13 

5.3 

1.45 

0.036 

0.276 

11-14 

4.4 

1.51 

0.043 

0.278 

*The  run  number  is  followed  by  a  coordination  number  which  is  used 
to  identify  data  as  recorded  on  magnetic  tape. 
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TABLE   2    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 

x  102 

2  X 

V 

12-1 

3.7 

2.49 

0.090 

0.346 

12-2 

5.6 

2.29 

0.063 

0.332 

12-3 

8.1 

2.13 

0.047 

0.250 

12-4 

10.9 

2.24 

0.060 

0.274 

12-5 

13.1 

2.13 

0.047 

0.218 

12-6 

15.5 

2.03 

0.053 

0.256 

12-7 

18.5 

2.03 

0.050 

0.256 

12-8 

19.2 

2.39 

0.048 

0.206 

12-9 

17.2 

2.08 

0.050 

0.252 

12-10 

14.1 

2.08 

0.048 

0.222 

12-11 

12.0 

2.18 

0.050 

0.246 

12-12 

9.3 

2.18 

0.057 

0.298 

12-13 

6.8 

2.18 

0.050 

0.268 

12-14 

4.8 

2.34 

0.077 

0.336 
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TABLE   2    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 
x   102 

A 

2  X 
V 

13-1 

3.6 

2.96 

0.120 

0.388 

13-2 

5.6 

2.65 

0.077 

0.376 

13-3 

8.8 

2.44 

0.053 

0.276 

13-4 

12.0 

2.65 

0.077 

0.332 

13-5 

15.4 

2.44 

0.053 

0.240 

13-6 

18.2 

2.39 

0.060 

0.290 

13-7 

19.1 

2.86 

0.055 

0.230 

13-8 

16.9 

2.34 

0.060 

0.282 

13-9 

13.8 

2.49 

0.057 

0.266 

13-10 

10.8 

2.65 

0.080 

0.342 

13-11 

9.4 

2.49 

0.063 

0.294 

13-12 

7.7 

2.44 

0.053 

0.262 

13-13 

6.6 

2.54 

0.060 

0.290 

13-14 

4.5 

2.70 

0.100 

0.375 
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TABLE   2    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 
(VRMS) 
x  102 

2  2 
V 

14-1 

3.6 

0.73 

0.015 

0.204 

14-2 

5.8 

0.73 

0.017 

0.204 

14-3 

8.3 

0.73 

0.016 

0.181 

14-4 

11.3 

0.68 

0.022 

0.219 

14-5 

14.1 

0.68 

0.016 

0.142 

14-6 

16.6 

0.68 

0.017 

0.156 

14-7 

20.5 

0.68 

0.015 

0.134 

14-8 

18.3 

0.68 

0.016 

0.148 

14-9 

15.2 

0.73 

0.017 

0.142 

14-10 

12.7 

0.73 

0.018 

0.158 

14-11 

9.7 

0.73 

0.017 

0.188 

14-12 

7.2 

0.73 

0.016 

0.188 

14-13 

6.2 

0.73 

0.017 

0.204 

14-14 

5.3 

0.73 

0.016 

0.204 

14-15 

4.2 

0.73 

0.016 

0.220 
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TABLE   2    (Cont'd) 


A 

Run 

f 

q 

eA 

2  y- 

V 

(CPS) 

(psf) 

(VRMS) 
x   102 

15-1 

3.6 

3.28 

0.140 

0.428 

15-2 

6.2 

2.91 

0.090 

0.378 

15-3 

9.9 

2.70 

0.598 

0.246 

15-4 

12.5 

2.86 

0.080 

0.316 

15-5 

15.1 

2.76 

0.060 

0.256 

15-6 

17.2 

2.54 

0.070 

0.294 

15-7 

20.8 

3.28 

0.058 

0.232 

15-8 

18.5 

2.60 

0.070 

0.292 

15-9 

16.2 

2.60 

0.066 

0.264 

15-10 

13.8 

2.76 

0.063 

0.272 

15-11 

11.3 

2.96 

0.090 

0.336 

15-12 

8.8 

2.65 

0.053 

0.258 

15-13 

7.5 

2.76 

0.066 

0.298 

15-14 

4.9 

3.01 

0.120 

0.390 
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TABLE  3 
REDUCED  DATA:  LARGE  VANES 


/\ 

f 

q 

e. 

2  2 

Run* 

A 

V 

(CPS) 

(psf) 

(VRMS) 
x   102 

5-1 

3.7 

2.29 

0.150 

1.084 

5-2 

5.6 

2.08 

0.120 

0.978 

5-3 

7.8 

1.40 

0.080 

0.922 

5-4 

10.0 

1.71 

0.110 

0.894 

5-5 

12.5 

1.56 

0.100 

0.876 

5-6 

14.3 

1.45 

0.100 

0.776 

5-7 

17.0 

1.04 

0.120 

1.072 

5-8 

19.3 

2.08 

0.093 

0.704 

5-9 

16.0 

1.04 

0.120 

0.960 

5-10 

13.0 

1.56 

0.090 

0.876 

5-11 

9.3 

1.45 

0.079 

0.776 

5-12 

6.6 

1.66 

0.170 

0.910 

5-13 

4.7 

1.77 

0.130 

1.118 

5-14 

3.6 

2.24 

0.150 

1.108 

*The  run  number  is  followed  by  a  coordination  number  which  is  used 
to  identify  data  as  recorded  on  magnetic  tape. 
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TABLE  3    (Cont'd) 


£i 

A 

o    V 

Run 

r 
(CPS) 

q 

(psf) 

eA 
(VRMS) 
x   102 

2  — 
V 

6-1 

3.6 

1.97 

0.145 

1.040 

6-2 

5.2 

1.56 

0.120 

1.186 

6-3 

7.5 

1.40 

0.077 

0.870 

6-4 

10.0 

1.51 

0.100 

0.842 

6-5 

12.4 

1.48 

0.101 

0.908 

6-6 

14.6 

1.30 

0.110 

0.906 

6-7 

16.8 

1.04 

0.120 

1.138 

6-8 

19.4 

1.77 

0.090 

0.584 

6-9 

15.2 

1.20 

0.108 

1.004 

6-10 

13.6 

1.40 

0.090 

0.816 

6-11 

11.3 

1.56 

0.110 

0.928 

6-12 

8.3 

1.25 

0.082 

0.982 

6-13 

6.5 

1.56 

0.090 

1.032 

6-14 

4.4 

1.66 

0.120 

1.200 
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TABLE  3    (Cont'd) 


K 

Run 

f 
(CPS) 

q 

(psf) 

eA 
(VRMS) 
x  102 

2  Z 
V 

0.130 

7-1 

3.7 

1.92 

0.924 

7-2 

6.4 

1.56 

0.100 

1.076 

7-3 

8.6 

1.25 

0.077 

0.850 

7-4 

10.3 

1.30 

0.083 

0.784 

7-5 

13.0 

1.14 

0.097 

0.766 

7-6 

15.6 

1.20 

0.097 

0.766 

7-7 

17.7 

1.04 

0.110 

0.986 

7-8 

21.0 

1.14 

0.083 

0.626 

7-9 

16.2 

1.04 

0.110 

0.876 

7-10 

14.2 

1.35 

0.083 

0.672 

7-11 

11.5 

1.45 

0.100 

0.832 

7-12 

9.3 

1.14 

0.081 

0.784 

7-13 

7.5 

1.35 

0.077 

0.816 

7-14 

5.1 

1.45 

0.110 

1.062 
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TABLE   3    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 
x  102 

A 
2  V 

V 

8-1 

4.0 

1.61 

0.109 

0.864 

8-2 

6.5 

1.40 

0.087 

0.922 

8-3 

8.4 

1.19 

0.066 

0.708 

8-4 

10.6 

1.14 

0.083 

0.724 

8-5 

13.2 

1.25 

0.083 

0.652 

8-6 

16.0 

1.14 

0.083 

0.638 

8-7 

21.0 

1.30 

0.077 

0.480 

8-8 

18.3 

1.46 

0.110 

0.874 

8-9 

14.1 

1.14 

0.077 

0.596 

8-10 

12.3 

1.30 

0.090 

0.720 

8-11 

9.5 

1.14 

0.073 

0.682 

8-12 

7.5 

1.25 

0.070 

0.734 

8-13 

5.2 

1.30 

0.090 

0.920 
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TABLE  3    (Cont'd) 


Run 

f 
(CPS) 

q 

(psf) 

eA 
(VRMS) 
x  102 

2  2 
V 

9-1 

3   6 

0   68 

fl   047 

0.670 
0.618 

9-2 

6.1 

0.62 

0.050 

9-3 

9.6 

0.57 

0.063 
0.043 

0.702 

9-4 

11.8 

0.62 

0.412 

9-5 

14.7 

0.52 

0.066 

0.622 

9-6 

18.1 

0.47 

0.073 

0.674 

9-7 

20.9 

0.52 

0.044 

0.424 

9-8 

16.7 

0.47 

0.059 

0.626 

9-9 

13.4 

0.62 

0.039 

0.388 

9-10 

10.6 

0.62 

0.060 

0.644 

9-11 

8.5 

0.57 

0.043 

0.540 

9-12 

7.3 

0.62 

0.047 

0.594 

9-13 

5.5 

0.68 

0.050 

0.670 

9-14 

4.4 

0.68 

0.043 

0.670 
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TABLE  4 

Experimental  Phase  Difference  Between 

Anemometer  and  Hot  Wire 

Small  Vanes 


q  =  4. 

.19  psf 

q  =  2.34  psf 

Run 

f 

<*> 

Run 

f 

4> 

(CPS) 

(degrees) 

(CPS) 

(degrees) 

1-1 

3.6 

7 

2-1 

3.6 

*  7 

1-2 

5.0 

12 

2-2 

6.0 

13 

1-3 

7.7 

15 

2-3 

8.8 

21 

1-4 

10.0 

20 

2-4 

11.2 

12 

1-5 

13.0 

17 

2-5 

13.3 

22 

1-6 

16.3 

22 

2-6 

15.0 

24 

1-7 

19.2 

27 

2-7 

16.9 

24 

1-8 

17.8 

22 

2-8 

19.1 

27 

1-9 

15.2 

23 

2-9 

18.1 

29 

1-10 

11.6 

18 

2-10 

15.9 

28 

1-11 

9.2 

16 

2-11 

14.1 

22 

1-12 

6.4 

12 

2-12 

10.1 

12 

1-13 

4.2 

9 

2-13 

7.0 

13 

1-14 

3.6 

6 

2-14 

5.0 

12 

50 


TABLE  4    (Cont'd) 


q  -  7, 

.55  psf 

q  =  5.90  psf 

f 

0 

Run 

Run 

f 

0 

(CPS) 

(degrees) 

(CPS) 

(degrees) 

3.6 

11 

3-1 

4-1 

3.6 

3 

3-2 

5.8 

15 

4-2 

5.1 

7 

3-3 

8.6 

21 

4-3 

7.9 

8 

3-4 

11.1 

29 

4-4 

10.1 

10 

3-5 

13.6 

32 

4-5 

12.0 

18 

3-6 

17.4 

36 

4-6 

14.0 

17 

3-7 

19.2 

32 

4-7 

17.0 

21 

3-8 

15.7 

37 

4-8 

19.3 

15 

3-9 

12.4 

32 

4-9 

15.6 

21 

3-10 

9.7 

23 

4-10 

13.1 

19 

3-11 

7.0 

20 

4-11 

9.1 

9 

3-12 

4.8 

15 

4-12 

6.6 

8 

3-13 

4.2 

13 

4-13 

4.6 

4 

4-14 

4.0 

4 
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TABLE  4    (Cont'd) 
q   =   1.00  psf 


Run 

f 

0 

(CPS) 

(degrees) 

10-1 

3.6 

23 

10-2 

5.5 

29 

10-3 

7.4 

33 

10-4 

9.7 

29 

10-5 

12.3 

36 

10-6 

14.5 

38 

10-7 

16.7 

49 

10-8 

19.2 

56 

10-9 

17.6 

48 

10-10 

13.3 

37 

10-11 

10.9 

34 

10-12 

8.4 

30 

10-13 

6.3 

29 

10-14 

4.3 

26 
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TABLE  5 
Experimental  Phase  Difference  Between 
Anemometer  and  Hot  Wire 
Medium  Vanes 


q  =  1.43  psf 


2.23  psf 


Run 

f 

0 

Run 

f 

0 

(CPS) 

(degrees) 

(CPS) 

(degrees) 

11-1 

3.7 

1 

12-1 

3.7 

1 

11-2 

6.4 

3 

12-2 

5.6 

2 

11-3 

8.6 

3 

12-3 

8.1 

2 

11-4 

11.2 

1 

12-4 

10.9 

0 

11-5 

13.1 

11 

12-5 

13.1 

5 

11-6 

15.5 

46 

12-6 

15.5 

40 

11-7 

18.2 

30 

12-7 

18.5 

23 

11-8 

19.3 

27 

12-8 

19.2 

24 

11-9 

16.8 

30 

12-9 

17.2 

22 

11-10 

14.2 

37 

12-10 

14.1 

23 

11-11 

9.7 

1 

12-11 

12.0 

0 

11-12 

7.7 

0 

12-12 

9.3 

0 

11-13 

5.3 
4.4 

3 
0 

12-13 

6.8 

0 

11-14 

12-14 

4.8 

0 
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TABLE  5    (Cont'd) 


q  =  2.56  psf 

q  m   0.71  psf 

Run 

f 
(CPS) 

0 

(degrees) 

Run 

f 
(CPS) 

0 

(degrees) 

1 

14-1 

3.6 

13-1 

3.6 

16 

13-2 

1  O   O 

5.6 

O   O 

3 
0 

14-2 
14-3 

5.8 
8.3 

15 

7 

13-3 

8.  8 

13-4 

12.0 

0 

14-4 

11.3 

8 

13-5 

15.4 

18 

14-,5 

14.1 

20 

13-6 

18.2 

28 

14.6 

16.6 

42 

13-7 

19.1 

22 

14-7 

20.5 

34 

13-8 

16.9 

36 

14.8 

18.3 

27 

13-9 

13.8 

0 

14-9 

15.2 

44 

13-10 

10.5 

0 

14-10 

12.7 

2 

13-11 

9.4 

0   . 

14rll 

9.7 

21 

13-12 

7.7 

1 

14-12 

7.2 

12 

13-13 

6.6 

0 

14-13 

6.2 

13 

13-14 

4.5 

1 

14-14 

5.3 

15 

14-15 

4.2 

14 
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TABLE  5    (Cont'd) 
q   =   2.85   psf 


Run 

f 
(CPS) 

(degrees) 

15-1 

3.6 

1 

15-2 

6.2 

2 

15-3 

9.9 

0 

15-4 

12.5 

0 

15-5 

15.1 

4 

15-6 

17.2 

40 

15-7 

20.8 

20 

15-8 

18.5 

27 

15-9 

16.2 

33 

15t10 

13.8 

0 

15-11 

11.3 

0 

15-12 

8.8 

1 

15-13 

7.5 

0 

15-14 

4.9 

0 
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Note: 


TABLE  6 
COMPARISON  OF  FOURIER  COMPONENTS 


Run 

S 
(x  102) 

Signal 

Fundamental 

0i 

2nd 
A2/A 

Harmonic 
1  </>l 

3rd  Harmonic 
A3^A1  02> 

(degrees) 

<%) 

<d, 

egrees) 

(%)  (degrees) 

1-1 

3.47 

eH 

— 

19.7 

- 

4.4 

- 

eA 

-2.0 

8.5 

-13 

10.7 

47 

1-2 

4.83 

eH 

- 

28.7 

- 

3.4 

- 

eA 

7.3 

30.8 

56 

2.8 

23 

1-2* 

4.83 

eH 

- 

30.5 

- 

6.0 

- 

eA 

7.5 

24.6 

20 

9.4 

-69 

1-3 

7.45 

eH 

- 

27.4 

- 

2.1 

- 

eA 

14.0 

8.1 

5 

6.4 

139 

1-4 

9.69 

eH 

- 

10.0 

- 

2.8 

- 

eA 

7.0 

9.0 

-104 

1.6 

162 

1.  (J>,      n  is  phase  difference  between  hot  wire  and  anemometer 
signal  for  component  in  subscripts  (positive  for  hot 
wire  leading  anemometer) . 


2     (  )  indicates  amplitude  ratio  of  components  in  subscripts 

A  for  specified  signal,  i.e.,  eH  or  e.. 

(  )                             n     a 

*  indicates  different  waveform  sample  in  run. 
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TABLE  6  (Cont'd) 


Run       S       Signal    Fundamental   2nd  Harmonic   3rd  Harmonic 
(x  102)  0.        A2/Ax  ^  A3/A]_   ^3 


I 


(degrees)  (%)    (degrees) 


(degrees) 


2-1  4.69 


2-14  6.15 


2-13  8.51 


2-12        12.25 


2-4  14.55 


2-4*        14.55 


2-11  18.4 


CH 
eA 
eH 
eA 
eH 
eA 
eH 
eA 
eH 
eA 

eH 
eA 
eH 
eA 


14 


28 


0.2 


12 


20.9 

— 

1.1 

— 

16.6 

10 

8.3 

56 

27.2 

- 

4.8 

- 

21.3 

-13 

9.4 

46 

23.4 

- 

2.8 

- 

8.4 

50 

9.4 

5 

26.5 

- 

4.4 

- 

8.8 

57 

7.8 

-16 

16.7 

- 

2.8 

- 

15.7 

-27 

1.6 

45 

16.2 

- 

5.1 

- 

16.3 

25 

1.1 

-25 

24.5 

- 

4.3 

- 

35.3 

1 

2.3 

-14 
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TABLE  7 
REDUCED  DATA:  SMALL  VANES 


A 

Run 

f 

q 

eA 

2  2 

V 

(CPS) 

(psf) 

(VRMS) 
x  102 

16-1 

3.6 

4.32 

0.0310 

0.0723 

16-2 

5.4 

4.26 

0.0240 

0.0562 

16-3 

7.0 

4.26 

0.0200 

0.0474 

16-4 

8.4 

4.26 

0.0170 

6. 0444 

16-5 

9.4 

4.26 

0.0180 

0.0526 

16-6 

10.8 

4.26 

0.0200 

0.0621 

16-7 

11.9 

4.26 

0.0140 

0.0503 

16-8 

13.6 

4.26 

0.0104 

0.0416 

16-9 

14.2 

4.26 

0.0104 

0.0397 

16-10 

15.6 

4.26 

0.0104 

0.0387 

16-11 

16.6 

4.26 

0.0096 

0.0387 

16-12 

17.7 

4.26 

0.0097 

0.0399 

16-13 

19.1 

4.26 

0.0095 

0 . 0444 

16-14 

18.3 

4.26 

0.0096 

0.0411 

16-15 

16.3 

4.26 

0.0099 

0.0387 

16-16 

13.7 

4.26 

0.0104 

0.0406 

16-17 

11.7 

4.26 

0.0148 

0.0511 

16-18 

9.3 

4.26 

0.0191 

0.0515 

16-19 

7.4 

4.26 

0.0180 

0.0454 

16-20 

6.1 

4.26 

0.0212 

0.0507 
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TABLE   7    (Cont'd) 


f 
(CPS) 

q 

(psf) 

A 

Run 

eA 

(VRMS) 
x  102 

2  X 
V 

3.6 

6.03 

17-1 

0.0524 

0.0793 

17-2 

5.3 

5.98 

0.0392 

0.0609 

17-3 

6.8 

5.98 

0.0320 

0.0507 

17-4 

8.1 

5.98 

0.0279 

0.0473 

17-5 

9.5 

5.98 

0.0341 

0.0588 

17-6 

10.4 

5.98 

0.0382 

0.0684 

17-7 

12.0 

5.98 

0.0247 

0.0527 

17-8 

13.1 

5.98 

0.0183 

0.0456 

17-9 

14.2 

5.98 

0.0150 

0.0429 

17-10 

15.7 

5.92 

0.0139 

0.0411 

17-11 

18.5 

5.98 

0.0139 

0.0470 

17-12 

19.1 

5.98 

0.0139 

0.0489 

17-13 

18.0 

5.98 

0.0139 

0.0454 

17-14 

16.6 

5.92 

0.0139 

0.0428 

17-15 

13.2 

5.98 

0.0183 

0.0451 

17-16 

11.0 

5.98 

0.0341 

0.0647 

17-17 

9.0 

5.98 

0.0299 

0.0533 

17-18 

7.0 

5.98 

0.0310 

0.0505 
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TABLE   7    (Cont'd) 


Run 

f 

q 

eA 

A 

2  2 

V 

(CPS) 

(psf) 

(VRMS) 
x  102 

18-1 

3.6 

7.60 

0.0784 

0.0836 

18-2 

6.4 

7.53 

0.0469 

0.0562 

18-3 

6.8 

7.53 

0.0428 

0.0540 

18-4 

8.1 

7.51 

0.0398 

0.0499 

18-5 

8.4 

7.51 

0.0377 

0:0507 

18-6 

9.3 

7.53 

0.0449 

0.0584 

18-7 

10.7 

7.53 

0.0531 

0.0709 

18-8 

11.4 

7.48 

0.0459 

0.0640 

18-9 

12.2 

7.48 

0.0346 

0.0546 

18-10 

12.8 

7.48 

0.0314 

0.0499 

18-11 

13.1 

7.48 

0.0293 

0.0491 

18-12 

14.0 

7.53 

0.0251 

0.0454 

18-13 

15.7 

7.48 

0.0173 

0 . 0449 

18-14 

17.2 

7.48 

0.0162 

0.0465 

18-15 

19.1 

7.60 

0.0173 

0.0554 

18-16 

18.3 

7.53 

0.0162 

0.0502 

18-17 

15.5 

7.48 

0.0173 

0.0445 

18-18 

13.7 

7.48 

0.0251 

0.0463 

18-19 

11.1 

7.48 

0.0490 

0.0669 

18-20 

8.3 

7.48 

0.0387 

0.0497 
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TABLE  8 
REDUCED  DATA:  MEDIUM  VANES 


Run 

f 
(CPS) 

q 

(psf) 

eA 

(VRMS) 
x  102 

2  X 
V 

19-1 

3.6 

2.34 

0.084 

0.370 

19-2 

5.5 

2.18 

0.058 

0.326 

19-3 

7.1 

2.08 

0.043 

0.238 

19-4 

8.6 

2.05 

0.041 

0.230 

19-5 

9.6 

2.08 

0.052 

0.266 

19-6 

11.3 

2.13 

0.055 

0.290 

19-7 

12.7 

2.08 

0.044 

0.222 

19-8 

14.2 

2.03 

0.046 

0.216 

19-9 

15.8 

1.98 

0.050 

0.246 

19-10 

17.6 

2.08 

0.046 

0.236 

19-11 

19.1 

2.34 

0.040 

0.206 

19-12 

18.0 

2.13 

0.046 

0.232 

19-13 

17.2 

2.03 

0.047 

0.238 

19-14 

15.2 

1.98 

0.049 

0.236 

19-15 

13.8 

2.08 

0.044 

0.214 

19-16 

13.0 

2.08 

0.043 

0.216 

19-17 

11.9 

2.08 

0.048 

0.248 

19-18 

10.6 

2.18 

0.060 

0.312 

19-19 

7.8 

2.03 

0.040 

0.226 
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TABLE  8    (Cont'd) 


y\ 

Run 

f 

q 

eA 

2  X 
V 

(CPS) 

(psf) 

(VRMS) 
x  102 

20-1 

3.6 

3.33 

0.138 

0.440 

20-2 

6.9 

2.83 

0.064 

0.270 

20-3 

7.7 

2.70 

0.056 

0.248 

20-4 

9.0 

2.81 

0.057 

0.240 

20-5 

10.6 

3.02 

0.091 

0.352 

20-6 

12.5 

2.86 

0.062 

0.258 

20-7 

13.5 

2.83 

0.057 

0.236 

20-8 

15.0 

2.70 

0.065 

0.254 

20-9 

16.4 

2.70 

0.066 

0.276 

20-10 

18.4 

3.17 

0.061 

0.258 

20-11 

19.0 

3.54 

0.057 

0.238 

20-12 

17.7 

2.91 

0.064 

0.270 

20-13 

14.1 

2.81 

0.059 

0.234 

20-14 

11.6 

2.86 

0.072 

0.300 

20-15 

8.0 

2.70 

0.054 

0.238 
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O-  LATERAL    AXIS 
A~  DRAG   AXI  S 
BRIDGE    EXCITATION  =  10.0  VOLTS(D.C) 


FIG.  2:    ANEMOMETER    STATIC    CALIBRATION 
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FIG.  3:  ANEMOMETER     INERTIAL    RESPONSE 
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FIG.  4:    SHUTTER   VANE    BLOCKAGE 
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FIG. 7:    SENSOR    INSTALLATION    DETAILS 
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FIG. 8:  ANEMOMETER    DRAG    vs.  DYNAMIC 
PRESSURE   (W.C.R.  TUNNEL) 
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FIG.  9:  MEAN  ANEMOMETER     DRAG    vs.  MEAN 

DYNAMIC     PRESSURE  (NONSTEADY   TUNNEL) 
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FIG.  10:  UNSTEADY    DRAG  vs.  FREQUENCY 
(SMALL  VANES) 
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FIG.  11:  UNSTEADY   DRAG  vs.  FREQUENCY 
(MEDIUM    VANES) 
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FIG. 12:  UNSTEADY  DRAG  vs.  FREQUENCY 
(LARGE  VANES) 
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** 


FIG.  13:  ANEMOMETER    PHASE    LAG  (SMALL  VANES) 
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FIG.14:ANEMOMETER    PHASE   LAG  (MEDIUM  VANES) 
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FIG.27:DRAG  RATIO   VARIATION  (Medium  Vanes,  q =2. 8psf) 
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FIG.28:DRAG  RATIO  VARIATION ( Large  Vanes,q=0.6psf) 
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APPENDIX  A 
Accelerometer  Calibration: 

Calibration  of  the  Statham  accelerometer  was  obtained  by  measur- 
ing amplitude  and  frequency  of  motion,  plus  phase  shift,  between 
motion  and  accelerometer  output  when  the  accelerometer  was  mounted 
on  the  table  of  a  Calidyne  shaker  unit,  model  A88. 

Amplitude  of  motion  was  obtained  using  a  Bently  distance  detector 
unit  which  is  effectively  a  reluctance  gage  that  senses  the  change  in 
distance  of  the  air  gap  between  a  flat,  non-magnetic  conducting  mater- 
ial and  the  sensor  coil  face.   For  the  calibration  range  considered, 
the  reluctance  gage  frequency  response  was  the  same  as  the  static 
sensitivity  with  a  distance  resolution  of  +  20  micro-inches.   Frequency 
was  sensed  to  an  accuracy  of  +  0.1  cycles  per  second  (cps)  using  a 
conventional  electronic  counter  unit,  Hewlett-Packard,  model  522. 
Phase  angle  difference  was  sensed  by  an  AD-YU,  type  405  phase  meter 
which  had  an  accuracy  of  +  1.0  degree. 

The  Statham  accelerometer  was  operated  at  a  10.0  volt  (D.C.) 
bridge  excitation  level  and  had  an  approximate  static  sensitivity  of 
4.0  millivolts  per  "g"  where  the  gravitational  unit  "g"  corresponds 
to  a  standard  value  of  32.2  feet  per  second^.   The  dynamic  response 
of  the  accelerometer  is  shown  on  Figure  35  including  both  gain,  GC(6J  ), 
and  phase  lag  characteristic  A(fe,.   The  acceleration  response  was 
checked  for  linearity  at  constant  frequency  up  to  two  g's.   The  un- 
damped resonant  frequency,  fjr,  was  evaluated  as  being  218  cps  using 
the  criteria  of  ninety  degrees  phase  lag.   The  accelerometer  was 
estimated  to  have  a  damping  ratio  of  0.558  relative  to  critical  damp- 
ing, and  comparison  is  shown  on  the  frequency  response  plot. 
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The  calibration  may  be  considered  as  related  to  true  acceleration 
by  the  relationship: 
If  true  acceleration  is 

ii(t)  -  Asm  cjt 
Then  the  accelerometer  reading  will  be  (A. 01) 

ec(t)   =  AGcSm(^f-A<£c) 

Where  Gp  and  A<pc   are  obtained  from  Figure  32. 
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APPENDIX  B 
Anemometer  Corrections  For  Tunnel  Floor  Accelerations 

It  is  the  purpose  of  this  section  to  illustrate  a  method  for  cor- 
recting the  drag  sphere  anemometer  balance-voltage  reading  for  the 
influence  of  tunnel  floor  accelerations  in  the  horizontal  axial  di- 
rection.  A  sketch  of  the  environment  is  indicated  below: 
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Nonsteady  Wind  Tunnel  (Not  to  Scale) 
At  the  reference  velocity  station,  the  hot  wire  provides  an  indi- 
cation of  the  air  velocity.   This  station  will  be  used  as  a  source  for 
the  origin  of  time. 

V+  v(-t)   =  V  [  I  +  <T  Sin  cjt]  (B.oi) 

where  0"  represents  the  magnitude  of  velocity  perturbation  relative  to 
the  steady  velocity  value.   The  floor  will  experience  an  acceleration 
due  to  the  aerodynamic  loading  induced  by  the  operation  of  the  rotating 
shutter  mechanism,  and  this  will  occur  principally  in  the  horizontal 
axial  direction.   For  purposes  of  this  discussion,  let  us  assume  that 
the  motion  is  harmonic  in  character  at  the  dominant  forcing  frequency. 
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The  floor  acceleration  can  be  expressed  as: 

X(t)    =  C\   sin  (djt-  <t>)  (b.02) 

where  it  is  recognized  that  the  quantities  Ci  ,  and  <p , ,  are  not 
measurable  directly. 

The  Statham  accelerometer,  which  is  firmly  attached  to  the  floor 
beneath  the  drag  sphere  anemometer  installation,  was  aligned  to  re- 
spond to  accelerations  in  the  horizontal  axial  direction.   Its  elect- 
rical signal  output  due  to  the  acceleration  environment  of  equation 
B.02  will  be: 

ec(t)  ■=  C2  Sin(ojt-4>2)  (B.03) 

where  the  quantities  C£  and  <p ^   are  directly  measurable.   From  labora- 
tory shaker  table  calibrations,  it  is  known  that: 
Cz   ~     C,   Cac(cj)    ,  and 

A<£c  =  \  -  4 

Also  we  know  from  a  laboratory  shaker  table  calibration  of  the 
drag  sphere  anemometer  that 

AeA(t)  =  C3  sin(cot-<*>3) 

=   C3  Sinj/tot-4>)-A<k] 

where       C3  =    C.GJoj)    f      and      A (£A  =  <£3  -  </> 
Combining   equations   B.02   through  B.04  allows  us   to   estimate   an   inertial 
correction   to   the  drag   sphere  balance  reading.      Note   that   q) 1    is   can- 
celled out   and   is   thus  never  required. 

AeA(t)  =q^|^Sin[(^t-4>)-A<*>  +  AC*,]       (B.05) 

The  observed  drag  sphere  anemometer  balance  output  is  directly 
measured  and  may  be  expressed  as: 

e      =  Bsin(cjt-d)  (B.06) 
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and  it  is  desired  to  obtain  a  corrected  reading  in  the  form: 
eAcoK=  BCOKSin(^-<r) 

Solving  equation  B„07  in  terms  of  quantities  involved  in  equations 
B.05  and  B.06,  one  obtains: 


(B.07) 


B 


^2 


COR* 


nfc 


and 


V  =  arcfan 


COS 


Q    G 


*-^%COS(+2  +  A+A-H) 


(B.C8) 


It  should  be  noted  that  the  phase  relationship  between  the  un- 
steady velocity  and  drag  (hot  wire  and  anemometer)  is  the  correct 
relationship  between  the  velocity  at  the  hot  wire  and  the  drag  at  the 
anemometer.   However  it  is  not  the  relationship  between  the  velocity 
and  drag  at  the  same  point.   Due  to  the  longitudinal  separation  of 
hot  wire  and  anemometer  the  pressure  signals  producing  velocity  per- 
turbations reach  the  anemometer  prior  to  reaching  the  hot  wire.   Based 
on  the  speed  of  sound  in  air,  the  mean  flow  velocity  in  the  tunnel, 
the  frequency  of  oscillation,  and  the  distance  between  hot  wire  and 
anemometer,  the  maximum  phase  difference  between  velocity  signals  at 
the  hot  wire  and  anemometer,  for  the  range  of  test  conditions  used, 
would  be  7  degrees.   This  phase  difference,  due  to  the  above  factors, 
means  that  the  phase  angle  (or  time  lag)  between  a  velocity  perturbaticr. 
at  the  anemometer  and  the  anemometer  drag  response  is  greater  than 
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that  indicated  by  phase  angle  measurements  between  the  hot  wire  and 
the  anemometer.   Figure  36  shows  the  computed  phase  shift,  due  to  the 
axial  separation  of  the  hot  wire  and  anemometer  sensors, as  a  function 
of  frequency  for  four  representative  values  of  dynamic  pressure. 
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APPENDIX  C 
Hotwire  Response  to  Base  Acceleration 

An  elementary  beam  theory  analysis  has  been  made  for  the  canti- 
levered  rod  (supporting  the  hot  wire  probe  and  containing  electrical 
leads)  to  determine  the  velocity  of  the  probe  induced  by  inertial 
loads  during  tunnel  floor  vibrations  in  the  horizontal  plane.   The 
entire  length  of  the  rod  (shown  in  the  sketch  below)  was  assumed  to 
be  identical  to  the  hollow  stainless  steel  portion  in  specific  weight, 
cross  section,  and  elastic  properties. 
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SKETCH  OF  HOTWIRE  INSTALLATION 

The  length  of  the  rod  was  taken  to  be  6.5  inches  since  the  rod 
is  fixed  at  the  bottom  of  the  tunnel  floor  by  a  threaded  collar 
arrangement.   This  is  a  conservative  value  because  the  floor  pro- 
vides some  support  due  to  the  close  fit  between  the  rod  and  the  hole 
through  the  floor. 
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The  following  properties  and  values  were  chosen  or  computed 
for  the  analysis: 

Modulus  of  elasticity  E  =  29  x  106  lbs/in2 

Specific  Weight  t  =  0.284  lbs/in3 

Second  moment  of  the  area  with  respect  to  the  ZZ  axis; 


I   =  tt 

■"■zz 


■  JL   (d24"dl4) 
64 


Cross  sectional  area; 


Izz  =  1.895  x  10"4  iti4 


A  =  1   (d22-d!2) 


A  =  0.0396  in2 
I    =  6.5  in 
MA.'  =  1'A 

IU=   0.01125  lbs/in 
ia"=   0.05  aak 
MA  =  yOC*.'  +  ajul' 
MA.=   0.0118  lbs/in 
The  vertically  mounted  rod,  subjected  to  an  acceleration  of  one  "g" 
in  a  horizontal  direction,  is  equivalent  to  the  rod  mounted  horizon- 
tally, loaded  by  its  own  weight,  as  shown  in  the  following  sketch. 

.uniform  load  of  JAA.    lbs/in 
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Elementary  beam  theory  gives  the  deflection  of  the  free  end 

(y  =  A   )  as: 

.3 

x  =  W  A  (C.01) 

8  EI 

where,  W  =  wi   under  a  one  "g"  condition,  and  the  maximum  deflection 
is  then 

,3  ,_3 


x 


-  (0.0767)lbs  (6.5)J  in' 


°   8(2.9  x  106)  lbs/in2  (1.895  x  10"4)  in4 
xQ  =  0.00048  in 
If  the  load  were  assumed  to  be  applied  harmonically,  -i.e., 

W  =  WQ  sinc^t,  (C.02) 

then  the  deflections  of  the  free  end  can  be  approximately  expressed 
by 

x  =  xQ  sin  (jjt  (C.03) 

provided  the  frequencies  of  interest  are  small  compared  to  the  lowest 
natural  frequency  of  the  rod. 

The  natural  frequencies  of  a  continuous  uniform  beam  are  given  by 

where  for  a  cantilever,  Gi  =  1.875  is  the  lowest  eigenvalue. 
Mass  per  unit  length 

m  =  w/g 

m  -  3.06  x  10"5  lb  sec2/in2 
Substitution  of  values  into  equation  C.04  yields 


/29  x  106  lbs/in2  (1.895  x  104)  in4 
(3.06  x  10"5)lb  sec2/in^  (6.5)4  in4 

^l   =  1113  radians  per  sec 

or  f^  =   111   cycles  per  second 
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Since  the  maximum  frequency  of  interest  is  20  cycles  per  second, 
the  ratio  (f/f^)  is  sufficiently  small  for  an  order  of  magnitude 
analysis  based  on  (C.03): 

For  x  =  xQ  sin  cot,  the  velocity  of  the  probe  tip  is  x  =LJx     cos  ust, 

Thus  the  tip  velocity,  relative  to  the  base,  can  be  further  ex- 
pressed as  x  =  ncjx  cos  cot   where  n  is  the  number  of  "g's"  acceler- 
ation. 

Converting  units  from  radians  to  cycles  and  inches  to  feet  yields 
an  expression  for  the  maximum  tip  velocity  in  terms  of  the  test  varia- 
bles n  ("g"  factor)  and  f  (frequency);   x    =  6.36  nf  x  10"°  ft/sec 

During  the  nonsteady  tunnel  test  program  the  maximum  horizontal 
acceleration  recorded  was  one-fourth  "g"  (n  =  1/4),  while  maximum 
frequencies  were  20  cycles  per  second;  thus 

Xj^x  =  6.36  (1/4)  (20)  x  10~6  ft/sec 
xmax  =  31.8  x  10"6  ft/sec 

The  minimum  velocity  perturbations  which  occur  anywhere  in  the 
investigation  are  about  0.26  feet  per  second.   Therefore  x    repre- 
sents about  0.01  percent  and  is  negligible 


109 


APPENDIX  D 
Tunnel  Floor  Velocity  Effects  on  Hotwire  Signal 

The  following  analysis  was  used  to  determine  the  velocity  pro- 
duced at  the  hot  wire  as  a  result  of  tunnel  floor  velocity,  considering 
the  hot  wire  support  to  be  rigid. 

The  floor  motion  was  assumed  to  be  simple  harmonic,  and  expressed 
as 

XF  =  A  5/n  U)t  (D.01) 

from  which  the  floor  velocity  and  acceleration  are  given  by 

XF  =  CJA  cos  cut  <D-02> 


and 


respectively. 


XF=  -OJ   A  Sin  cjt  (D.03) 


2 
The  factor  do  A  can  be  determined  from  the  accelerometer  data 

for  each  run.   Then  using  the  oscillation  frequency  of  that  run  the 
maximum  displacement,  A,  can  be  determined,  and  the  product   A  is 
the  velocity  amplitude.   Due  to  the  rather  large  variations  in  accel- 
eration, and  the  fact  that  the  term  A  varies  with  frequency,  no  ex- 
pression can  be  obtained  which  describes  the  entire  test.   Instead 
the  quantity  A  must  be  computed  for  each  test  point. 
The  following  sample  calculation  indicates  the  method  used: 
Run  5-1: 

from  raw  data:    f  =  3.7  cycles  per  second 

ec  =  0.2  millivolts  (RMS) 
accelerometer  sensitivity:   kc  =  4.02  millivolts  per  "g" 
a/A  =  &c/k     =  0.0498  "g's"  (RMS) 

2, 


wA=  1.60  feet/sec2  (RMS) 
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Now  the  velocity  amplitude,  coA,    is  given  by 

X_  =coA  = 


_     6J2A  \.<oO 


F  GJ  £7ff 

Jip  ~  0.0689     H/5ec     (Rms) 

From  previous  data  reduction 

eH  =  0.  El  volts  (rms)  =>  v  =  24- -f-^ec  (rms) 

The  relative  magnitude  is  then 

X^  =  0.00287 

or  about  0.3  percent. 

The  maximum  value  of  x^/v   found  in  the  test  runs  was  of  the 
order  of  2  percent,  which  was  considered  to  be  a  negligible  influence. 
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APPENDIX  E 
Quasi-Steady  Theory: 

The  ability  of  the  drag  sphere  to  sense  dynamic  pressure  changes 
by  the  change  in  drag  force  is  clear  for  "steady"  velocity  changes, 
however,  the  ability  of  the  time-varying  drag  force  to  provide  a 
reasonable  representation  of  time-dependent  variation  of  dynamic 
pressure  cannot  be  assured  by  a  linearized  approach  unless  verified 
by  experiment.   The  quasi-steady  theory,  as  explained  here,  is  based 
upon  the  premise  that  the  drag  force  coefficient  remains  invariant 
and  that  the  instantaneous  value  of  drag  is  dependent  directly  upon 
the  instantaneous  value  of  dynamic  pressure.   Comparison  between  exper- 
iment and  quasi-steady  theory  is  the  purpose  of  this  research,  with 
the  goal  of  showing  both  the  drag  magnitude  and  phase  shift  as  a 
function  of  dimensionless  frequency.   Quasi-steady  theory  assumes  no 
phase  shift. 

If  the  velocity  is  assumed  to  be  a  harmonic  component  superim- 
posed upon  a  steady  term,  then  we  have  that: 

V+  v(t)    =  V(/  +  <7Sincjt)  (E.OD 

where  c  represents  the  amplitude  ratio  between  unsteady  and  steady 
velocity.   The  instantaneous  value  of  dynamic  pressure  can  then  be 
expressed  as: 

f4?V  [0  +  £)  +  2o-Sincjt-£-co«2cj*]   (EQ2) 

Applying  the  quasi-steady  assumption  to  drag  force  yields: 

D  +  d(t)--g-fCDAV2[0*f^+2<r  sin  o»t  -  f  cos  2  art]      (e.03) 
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It  will  be  noted  that  the  steady  drag  term  increases  slightly  over 
the  value  in  uniform  flow  due  to  the  presence  of  the  C  /2  term  in 
equation  E.03  .   As  explained  in  Section  3,  the  influence  of  the  term 
was  experimentally  verified  as  being  small,  which  is  reasonable  to 
expect  for  velocity  perturbations  up  to  about  15  percent  of  the  steady 
velocity;  i.e.,   |  v(t)|  <   0.15  V. 

From  equation  E.03,  we  can  readily  discern  that  the  ratio  of 
unsteady  drag  to  steady  drag  by  quasi-steady  theory  is  given  by: 

^  =  2o-smcut--|-2  cos  2cot  (e.04) 

Next  we  consider  equation  E,04  from  the  standpoint  of  mean  square 
values,  which  becomes  upon  applying  the  limit  operation  for  the  mean 
square  of  a  temporally  varying  term: 

-2.  =    Um    (±f(4o-2sin26jt  -2<r3sin^t  cos  2**  + 

Defining   d.  -*y  d    and  considering  equation  E.05, 

We  now  must  recognize  that  the  velocity  perturbation  O"  corresponds 
to  a  peak  amplitude  and  may  be  expressed  in  terms  of  the  RMS  of 
velocity  perturbation,  €     ,    by  0"  =  v£-  €    . 
Hence  we  obtain  that 

which  to  a  first  order  term  may  be  stated  as: 

4  =  2-$-  <E-06> 

Equation  E.06  is  a  statement  of  the  quasi-steady  drag  ratio  as  used 
in  Section  3  of  this  analysis. 
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APPENDIX  F 


Fourier  Analysis  of  Signal; 

A  simple  analysis  was  made  of  the  time  varying  electrical  signals 
from  the  hot  wire  and  drag  sphere  anemometer  in  order  to  evaluate  the 
quality  of  the  unsteady  velocity  perturbation  and  the  response  of  the 
anemometer  to  the  various  harmonics.   Strictly  speaking,  the  analysis 
using  a  Fourier  series  of  consecutive  harmonics  is  not  valid  for  this 
experimental  program,  since  it  is  possible  that  a  frequency  may  be 
present  in  the  signal  which  is  not  an  integer  multiple  of  £he  funda- 
mental, and  hence  a  Fourier  analysis  using  a  truncated  series  could 
give  a  false  indication  of  signal  content.   However,  the  Fourier 
approach  was  rationalized  in  this  instance  on  the  basis  that  the  time 
to  do  a  true  power  spectral  and  cross-spectral  analysis  was  not  con- 
sistent with  the  available  time,  and  at  least  an  approximate  indication 
of  frequency  content  had  to  be  obtained. 

Consider  a  periodic  function  of  period  2T  as  shown  in  the  sketch 
below: 
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The  Fourier  series  representation  for  the  cyclical  function  may 
be  expressed  as: 

f(x)  =  *°  +  2  (Ancos  *Z£  +  Bnsin  £2p)       <F.oi) 
2rn=l 

An  =  ^F  f  to)C05  ^^X   J   n«-Oyli^...      (F.02) 

r\  =  1/2,1,  •••     (f.03) 


where 

One  complete  cycle  of  the  signal  to  be  analyzed  was  obtained 
from  an  oscilloscope  photograph  (cf.,  Figure  30)  and  divided  into 
36  equal  parts,  such  that 

3&  A*X^  =  27r  radians   =  3CzO    degrees 

A^Cj  ~  10   degrees 

*i   =  101  degrees  (i«  O,  1,2, -",35) 

Now  express  period  as 

C T  =  3fc0  degrees 

Then  at  each  10  degree  station  (x.)  through  the  cycle,  the  amplitude, 

f(x^),  was  measured.   This  information  allowed  the  coefficients, 

A  and  B  ,  to  be  evaluated,  using  a  numerical  approximation  to  the 


integral  equations  F.02  and  F.03,  as 

i=o 
which,  with   x-  =  10  degrees,  is 

35 


An=  7^2  f(**)cM  nxi 


Similarly 


L=Q 


35 


Equation  (F.01)  may  be  alternately  expressed  as 

.  cQ 

-      hi 

where 


((x)--^+2Cncos(ur,t-4'r,) 


and 


4v,  =  arc  tan 


LA, 


(F.04) 

(F.05) 

(F.06) 

(F.07) 

(F.08) 
(F.09) 
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The  constant  term,  AQ,  is  a  function  of  the  axis  location  (from 
which  f(x-)  is  measured)  and  is  of  no  value  in  this  analysis.   Only 
three  components  were  computed,  i.e.,  the  fundamental,  second,  and 
third  harmonic.   After  Cn  was  computed  (n  =  1,2,3)  for  a  given  signal, 
the  magnitudes  of  the  second  and  third  harmonic  components,  relative 
to  the  fundamental,  was  given  by  the  ratios,  C2/C-^  and  Co/C^.   Compu- 
tation of  c£n  provided  phase  relationships  among  components  of  one 
signal,  and  also  phase  relationships  between  two  signals  processed 
relative  to  the  same  time  scale  and  origin. 
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